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This is a pedagogical introduction to Brownian motion on the occasion of the 100th anniversary of Einstein's 
1905 paper on the subject. After briefly reviewing Einstein's work in its contemporary context, we pursue 
some lines of further developments and applications in soft condensed matter and biology. Over the last 
century Brownian motion became promoted from an odd curiosity of marginal scientific interest to a guiding 
theme pervading all of the modern (live) sciences. 



1 Introduction: 100 years ago 

"If, in some cataclysm, all of scientific knowledge were to be destroyed, and only one sentence passed 
on to the next generation of creatures, what statement would contain the most information in the fewest 
words? I believe it is [. . . ] that all things are made of atoms — little particles that move around in perpetual 
motion, attracting each other when they are a little distance apart, but repelling upon being squeezed into 
one another." Some decisive theoretical tools to put the atomic hypothesis (here quoted from The Feynman 
Lectures [ 1 ] ) at scrutiny were provided in Einstein's paper "Uber die von der molekulartheoretischen Theo- 
rie der Warme geforderte Bewegung von in ruhenden Fliissigkeiten suspendierten Teilchen" [2] published 
in this journal 100 years ago, which starts with the claim "that according to the molecular-kinetic theory of 
heat, bodies of a microscopically visible size suspended in liquids must, as a result of thermal molecular 
motions, perform motions of such magnitudes that they can be easily observed with a microscope." It 
quickly became clear that Einstein's calculations could explain hitherto mysterious observations of dancing 
pollen granules and other small particles first systematically described by the Scottish botanist Robert 
Brown [3] and (independently) by his French colleague Adolphe Brongniart [4]. And a few years later, 
drawing upon painstaking observations of thousands of Brownian particles in his lab at the Sorbonne, Jean 
Perrin could ultimately end the debate about the reality of atoms [5]. Though dating back to the ancient 
Greeks, this debate had considerably intensified towards the turn of the century and prompted Planck's 
dictum that the acceptance of new scientific ideas proceeds via natural extinction rather than by persuasion 
of the opponents, which is sometimes referred to as Planck's principle among philosophers of science. 

It was soon realized that the laws of Brownian motion do not exclusively pertain to the realm of physics 
but likewise to chemistry, biology and even economics. After the dust had settled, this was acknowledged 
by a couple of Nobel prizes: in physics, Perrin was rewarded for proving the reality of atoms while the 
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chemists at the same time honored Svedberg for his dispersion studies and for the invention of the ul- 
tracentrifuge, and later Onsager for further growing the germ layed in Einstein's 1905 paper. Also, the 
economists eventually embraced the laws of Brownian motion as useful tool for the valuation of stock 
options in due time for the centenary of their derivation by Bachelier [6]. Since he (and other contributors) 
had unfortunately passed away by that time, the prize went to Merton and Scholes 1 . 

From the historical perspective, the theory brought forward in Ref. [2] (together with related work by 
Smoluchowski [8] and Langevin [9]) is therefore certainly not to be sniffed at by her two sisters, relativity 
and quantum mechanics. It took quite a different course, though. While the latter two were soon promoted 
to the holy grail of science and have never since ceased to attract the public, the former catched much less 
attention. This is all the more astonishing as the vision it promoted, namely that our macroscopic world 
including its living organisms resides on a jittering sea of fluctuations, is particularly remarkable for the 
broadness of its impact. During the last 100 years, the facets of this vision have continuously multiplied, 
thereby providing ample evidence for a "slower" [10], yet "unfinished" [11] revolution. Brownian mo- 
tion remains a leading paradigm in a development that seems directed oppositely to that of the quantum 
& relativity revolution — not deeper into the elements but towards the understanding of the principles 
underlying the emergence of collective properties in many-particle assemblies — and for which P. W. 
Anderson has coined the mantra "more is different" [14], meaning that "at each new level of complexity, 
entirely new properties appear, and the understanding of these behaviors requires research which I think is 
as fundamental in its nature as any other." Depending on the importance of thermal fluctuations one may 
roughly divide this research into "hard matter" ("/cbT — ► 0") and "soft matter" [15, 16] ("fc^T ~ l") 3 
physics. Underlying both is the discipline of statistical mechanics, which has not only become a major part 
of any physics curriculum, but which has over the years proven its vitality by invading numerous foreign 
disciplines such as biology, ecology, economics, traffic, sociology, game theory, etc. [19-24], and by iden- 
tifying new interdisciplinary research topics like criticality, structure formation and synergetics, quenched 
disorder, glassy dynamics, and many others, where it tries to bring to bear its concepts and methods. On the 
other hand, widespread hopes that on a more fundamental, sub-atomic scale fluctuations would eventually 
cease and give way to solidity and firm order did not quite materialize. Instead, the toolbox of statistical 
mechanics turned out to be instrumental also to various "fundamental" developments from quantum theory 
to cosmology (path integrals, stochastic quantization, space-time foam, . . . ). 

In the following, our scope will have to be much, much narrower, though, and we will try to restrict 
ourselves to tracing a selective set of implications and modern offsprings of Einstein's 1905 paper on 
Brownian motion [2], with the aim of illustrating variations of its basic themes in soft condensed matter 
and biological physics. Since, above all, Brownian motion is particularly strongly linked to biology, not 
only by its name; a link that we welcome as a guide (or excuse) for taking shortcuts through the bewildering 
jungle of applications. Soft matter and biological physics are both fields governed by meso-scale structures 
for which thermal fluctuations are all important, the relevant dynamics often bridging the gap between the 
microscopically fast and the macroscopically slow, as neatly summarized by the formulas 

fc B Tw4pN-nm, k B T/r) w w 4 ^mV 1 , (1) 

where ?y w denotes the viscosity of water. Conversely, these formulas immediately raise the (legitimate) 
expectation that Brownian motion is an ubiquitous phenomenon in these fields. 

More precisely, our outline is as follows. After a brief historical introduction to Einstein's paper [2] 
and closely related work in section 2, the remainder is organized along six traces. Our starting point 
is the observation that the derivation in Ref. [2] inspired and partly anticipated major developments in 

1 Commonly better known for their role in one of the world's largest financial collapses [7]. 

^ In view of beautiful recent progress in what one might call "condensed matter cosmology" (e. g. [12, 13]), this antithesis seems 

less compelling now than some years ago. 
3 Here, the characteristic intrinsic energy scale of the system is understood to provide the unit. In this sense, kgT f» 1 may for 

instance indicate that Brownian motion limits our auditory abilities [17], or the accuracy of a gravitational wave detector [18] 

or (in the form of Nyquist-Johnson noise) of any electrical measuring instrument. 
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non-equilibium statistical mechanics (section 3) as well as a powerful stochastic calculus (section 4) with 
innumerable applications, of which we only mention a few paradigmatic examples taken from biology in 
section 5. Precisely as proposed by Einstein in 1905 [2], Brownian particles have ever since served as 
versatile probes of their environment (section 6). Moreover, their stochastic fluctuations have provided 
valuable insights into their own (soft) internal degrees of freedom, which can give rise to most complex 
macroscopic material properties thus providing a unique potential for applications (section 7). Finally, 
Brownian fluctuations are both exploited and rectified in ingenious biological micro-machines (section 8). 

Needless to say that, even within this limited scope, we cannot be objective or complete in any sense. 
Original publications can by no means be cited systematically, and citations should rather be taken as our 
(certainly biased) personal suggestions for further reading. With these cautious remarks, we hope that the 
reader is prepared for a very subjective random walk lurching to and fro between physics and biology. 

2 Historical notes: from Brown to Perrin 

Brownian motion: a very brief history 

"The story of Brownian motion is one of confused experiment, heated philosophy, belated theory, and, 
finally, precise and decisive measurement" [25]. It is well documented in various publications by Perrin, 
Nye, Brush, and others, and summarized in several monographs on Brownian motion (e. g. [26-29]). Here, 
we content ourselves with a few remarks. First, since we expect 2005 to be a year with strong "Matthew 
effects" (in case of doubt, credit is ascribed to those who already got most [30]), it seems in order to start 
by a brief disclaimer: Brown was not the first observer of Brownian motion (similar observations had 
repeatedly been made after the first microscopes came up in the late 17th century); Perrin was not the first 
to show that kinetic theory is in accord with experiments, nor was the admirable diligence of his group 
necessary to make the atomic hypothesis acceptable to the majority of researchers; Einstein did of course 
not invent the diffusion equation — which he supposed to be generally known just as we would today — 
nor did he first suggest that Brownian motion is due to thermal motion of the solvent molecules or provide 
the first testable prediction of statistical mechanics or kinetic theory, etc. 

In fact, Brown was curious enough about the erratic dance under the microscope to spend some time 
in the late 1820's on systmatic investigations, trying to clarify its origin. Starting out with micrometer- 
sized granules from the pollen of Clarkia pulchella and later using several other substances, he convinced 
himself (though not everybody else) that internal animate as well as external perturbative causes and optical 
illusions could be ruled out. Nevertheless, the phenomenon hardly received the attention it deserved. It 
remained a curiosity during Brown's lifetime and almost until the end of the century, when Leon Gouy's 
detailed investigations of the phenomenon [31] led him to the conclusion that the ubiquitous perpetual 
motion is caused by thermal motion of the solvent molecules. (This was in 1888, when Einstein still 
went to primary school, and even at that time the idea was not really new 4 ) He later also noticed that this 
explanation seemed to be in contradiction to the second law of thermodynamics (a concern independently 
raised by Rontgen in a letter to Einstein in 1906) or could at least provide an ideal natural laboratory for 
direct experimental examination of certain apparent logical contradictions between the atomic hypothesis 
and thermodynamics. By that time diffusion was a well-established thermodynamic notion routed in 
the work "Uber Diffusion" (1855) by the German pathologist Adolf Fick, which starts with the words 
"Hydrodiffusion through membranes is not only one of the basic factors of organic life 5 , but also a most 
interesting physical process that deserves much more attention from physicists than it has got so far" [34]. 
But the relation between Brownian motion and diffusion — today often regarded as synonymous — was 



See e. g. C. Wiener, "Erklarung des atomistischen Wesens des tropfbar-fliissigen Korperzustandes und Bestatigung desselben 
durch die sogenannten Molecularbewegungen" [32]. 

Some advocates of the opposite view nowadays claim to have compelling evidence that "there is certainly no reason to suppose 
that classical diffusion theory, or any of its offspring, will play a significant role in our understanding of biology in the future". 
[33] 
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much less clear. This was the missing link provided by Einstein and Smoluchowski, who thereby pointed 
the way to quantitative hands-on experiences with a strange microcosmos awaiting discovery. Despite of 
work by Clausius, Maxwell, Boltzmann, Loschmidt and others on heat conduction and interdiffusion that 
had clearly demonstrated the usefulness of kinetic theory and even produced reasonably accurate estimates 
for the size of molecules by the late 19th century, there remained strong opposition against taking the 
underlying atomistic view too literal. Only Perrin's extensive studies, performed with Einstein's prediction 
at hand, could finally convince the last of the "Energetiker" around Ostwald that atoms were more than just 
a convenient analogy. To us, the great impact of Perrin's and Einstein's work on Brownian motion seems 
partly linked to the fact that, instead of only measuring another thermodynamic or hydrodynamic relation 
that could be rationalized by kinetic theory, a forthright manifestation of molecular chaos, namely the 
predicted "strange" equation of motion of individual diffusing colloids (on average the travelled distance 
grows with the square root of time), could directly be observed. The fractal trajectory implied by Einstein's 
result provided a convincing explanation why earlier attempts to measure the velocity of Brownian particles 
and compare it to the average thermal velocity predicted by kinetic theory (for free particles), were in vain: 
with the characteristic distance travelled by the particle growing as \/t, the velocity is proportional to 
Naturally, around 1900, a non-differentiable trajectory was hard to swallow, and Einstein hastened to add 
[35] that the molecular discreteness provides a lower cutoff to such unseemliness. Curiously, the same 
fractal trajectories were determined to reappear later in the disguise of the Heisenberg uncertainty relation, 
which indicated that they also prevail on a more fundamental level [36] (see section 4). 

Finally, it should be said that Einstein's interest in Brownian motion grew out of a profound interest 
in statistical fluctuations [35,37,38], which can be identified as a common denominator of his work in 
seemingly unrelated areas. In particular, fluctuations played a crucial role [39] in his persistent efforts to 
find a rational explanation for Planck's radiation law. Einstein realized that stochastic fluctuations rule on 
the quantum level (photon statistics) in classical small systems (Brownian motion) and near a critical point 
(critical opalescence), thereby anticipating at the beginning of the 20th century, what became vast fields of 
research towards its end. Without exaggerating, he may therefore be called the father of fluctuation theory 
[40]. 



Revision course: Einstein's 1905 paper on Brownian motion 

In brief, the argument in Ref. [2] goes as follows. If the atomic hypothesis holds, particles immersed in 
an ambient medium must undergo perpetual irregular motion due to chaotic collisions with its molecules. 
It was Einstein's insight (cited above) that micrometer-sized particles are both big enough to be seen 
themselves in the microscope and small enough that their Brownian motion is substantial, and, as a conse- 
quence, keeps them suspended indefinitely against gravity. The publication consists of two major parts. In 
the first part, Einstein argues against the view — still common at that time — that suspended particles may 
be dismissed in thermodynamics. He insists that N particles of arbitrary volume v suspended in a solvent 
at small number density n <C represent an ideal gas and therefore give rise to an osmotic pressure 
p = nksT (in modern notation), thereby excerting a finite, albeit small, thermodynamic force on their 
container. In presence of a constant volume force K acting on the particles one thus has the (isothermal) 
force balance n(r)K + Vp(r) = 0, or p(r) = feTii(r) oc e" K ' r / kBT in the steady state. Further, 
Einstein suggests to regard this stationary balance as a dynamic equilibrium between a diffusion current 
—DVn and a drift current nK/(, with ( and D the friction and diffusion coefficient of the Brownian 
particles in the solvent, respectively. Eliminating K and n between the balance of forces and that of drift 
and diffusion, he infers what is today known as Stokes-Einstein relation 



D = ksT/C 



(2) 
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So far, this is a relation between two coefficients appearing in two generalized hydrodynamic equations, 
namely the Stokes equation for incompressible viscous flow 6 , and the diffusion equation. This is not yet 
quite, what made it so remarkable. We must point out that the notion of Boltzmann's constant ks was 
not yet fully established at the time, and its apparence in Eq. (2), originally in the form of the gas constant 
divided by Avogadro's number, allowed the latter to be determined from measuring macroscopic quantities. 
Intriguingly, it thereby indicated a relation of the macroscopic kinetic coefficients D, ( to the microscopic, 
molecular world. 

This link, which was later formalized under the name of Green-Kubo relations, is established in the 
second part of Einstein's paper, which is quoted as a pedagogical introduction to stochastic processes 
in Gardiner's textbook [28]. It undertakes a probabilistic derivation of the diffusion equation based on 
a particularly simple model for the fluctuations of the Brownian particle, called a random walk. The 
essential idea is to express the particle concentration (in one space dimension) at position x and time t + t 
as a function of the concentration at time t, 



Due to the use of the auxiliary "jump propabilities" </v(0 in this Chapman-Kolmogorov equation, as 
it is called today, Einstein's derivation appears much more elegant and transparent than Smoluchowski's 
more mechanistic discussion of Brownian motion, because it circumvents the construction of an explicit 
microscopic dynamical theory of the process by some (implicit) technical assumptions on ip T (0- Taylor 
expanding n(x, t) with respect to the "small" quantities r and £ and matching the leading order terms on 
both sides, Einstein recovers the diffusion equation 



The expressions for D are understood to be independent of the considered time interval r as long as t is 
large compared to a microscopic collision time. In his derivation Einstein also invokes the limit "r — > 0", 
which later became the subject of debate (see below). The final relation between D and the mean square 
displacement (Sr 2 (t)) = ([r(t) — r(0)] 2 ) of a Brownian particle starting at t = at r = r(0) follows from 
applying the diffusion equation to each of the diluted particles separately. It is of interest as paradigm for 
a general rule: kinetic coefficients (here D, Q in macroscopic thermodynamic or hydrodynamic equations 
(here diffusion, Stokes equation) can entirely be expressed in terms of correlation functions of fluctuations 
of microscopic variables (see section 3). This is how the discrete chaotic processes on the microscale 
generate the smooth behavior on the macroscopic scale. 

3 Non-equilibrium statistical mechanics: from Stokes-Einstein relations 
to fluctuation theorems and effective temperatures 

Close to equilibrium: the fluctuation-dissipation connection 

The most celebrated result of Eintein's 1905 paper on Brownian motion [2] is the so-called Stokes-Einstein 
relation that links two hydrodynamic transport coefficients to each other and to thermal fluctuations, see 



This is nothing but the diffusion equation, Eq. (4), with D replaced by the kinematic viscosity, and rt(r) replaced by the solvent 
velocity field projected onto its "transverse" components (so that it is non-divergent). 




(3) 



d t n(r,t) = D\7 2 n(r,t) 



(4) 



along with a stochastic expression for the diffusion coefficient, 




(5) 
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Eqs. (2), (5). The second part of it, the relation between hydrodynamic transport coefficients and micro- 
scopic fluctuations, is today more commonly known as Green-Kubo relation, and is usually presented in 
the form [41] 



As Einstein himself noticed [35,37], the content of Eq. (6) is not restricted to Brownian motion but is 
straightforwardly translated to electrical circuits 7 , a case that was later further elaborated by Johnson and 
Nyquist [43,44]. Analogous relations have been derived for all kinds of transport coefficients (conductiv- 
ities, viscosities, etc.) relevant in different areas of condensed matter physics. They are of considerable 
practical interest, as they suggest a way to do non-invasive measurements of response coefficients (here 
Q. And they provide insight into the complicated behavior of strongly interacting many-particle systems. 
For example, the general form of Eq. (6) persists for dense suspensions of Brownian particles, where the 
gradient diffusion coefficient governing the diffusive spread of concenrtation variations via the diffusion 
equation becomes renormalized by the (osmotic) compressibility of the suspension compared to the bare 
diffusion coefficient of a single particle in the solvent. This knowledge allows one (among other things) 
to infer why the dynamics slows down near phase transitions, where susceptibilities (here the compress- 
ibility) diverge. Thus, near a critical point [45], fluctuations not only become long-ranged as foreseen by 
Einstein in his paper on critical opalescence [38], but they also become long-lived; though in a slightly 
more intricate way than suggested here [46,47]. 

It was Onsager who started the excavation of the treasure buried in the Stokes-Einstein relation. Follow- 
ing Einstein's later work [38] on fluctuations, he likened the relaxation of hydrodynamic variables in re- 
sponse to weak external perturbations to the decay of correlations between the corresponding microscopic 
variables: the forced deviation from equilibrium could as well have been a random fluctuation of the sys- 
tem. This regression principle allowed him to exploit exact microscopic symmetries, notably the invariance 
of the microscopic equations of motion under time reversal, to infer constraints on the kinetic coefficients: 
however difficult and incomprehensible e. g. the coefficient of thermodiffusion (Soret effect) 8 may seem, 
thanks to Onsager we can exactly relate it to a corresponding coefficient for "diffusive heating" (Dufour 
effect). Further generalizations of Eq. (2) like the mentioned Green-Kubo relations and the fluctuation- 
dissipation theorem followed. Around these cornerstones the field of linear non-equilibrium thermody- 
namics (also "irreversible thermodynamics", "generalized/molecular hydrodynamics") developed. The 
modern derivation [49-51] of all of the mentioned relations is based on linear-response theory (often 
complemented by projection operator techniques [52] to separate slow, hydrodynamic from fast, fluctuat- 
ing variables), which therefore can be regarded as the foundation of "close-to-equilibrium statistical me- 
chanics"; see Kubo's short pedagogical overview [41]. The basic assumption underlying linear-response 
theory is that the phase-space density may be linearized in a small perturbation added to the Hamilto- 
nian of the system. This has been the subject of some debate, reviewed in [53]: since single trajectories 
of nonlinear dynamical systems are known to exhibit sensitive dependence on perturbations, the required 
macroscopic linearity cannot be a property of single trajectories but rather has to rely on some ergodic 
property of the phase space dynamics as a whole [54]. Only for toy models such as the periodic Lorentz 
gas or Sinai billiard [55], the link between the transport coefficients and the underlying microscopic Hamil- 
tonian dynamics can be made fully explicit [53,56], while for the complex systems generally of interest 
in statistical mechanics this is out of reach. Already a rigorous microscopic derivation of Fourier's law 
for heat conduction, which is a diffusion equation with D replaced by the thermal conductivity, turns out 
to be quite an endeavor (most recently reviewed in Ref. [57]). Similarly, it was recently disputed (see 



Since he could not imagine other testable consequences, he dismissed any discussion of further realizations as "useless" [35] 
but maintained interest in the electrical analog [37], for which he even designed a "Maschinchen" (little machine) [42] to do 
experiments by himself. 

See e. g. Ref. [48] for a recent experimentally oriented review, and various contributions in Eur. Phys. J. E 15(1), 2004. 




oo 



dt (v(t) ■ v(0)} . 



(6) 
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[58]) whether extremely careful measurements [59] might be able to establish the role played by chaotic 
microscopic dynamics for Brownian motion. So the question, whether chaotic microscopic dynamics or 
the high dimensionality of the phase space is responsible for the irregularity of Brownian motion is, after a 
century, still open for investigation. 




Fig. 1 Probability distribution for the dissipative 
work W. The finite weight for negative VV implies 
that Einstein's theory of Gaussian fluctuations, if ap- 
plied to W, predicts transient violations of the second 
law. 



Fig. 2 Various scenarios for generalized fluctuation- 
dissipation relations represented in integrated form: 

X = tidrd(x(t))/dF(t') 



Far from equilibrium: "violating" basic laws and aging 

Far from equilibrium, the whole formalism of irreversible thermodynamics is a priori not generally appli- 
cable. However, it is a natural temptation to try to extend the approved formalism to classes of problems 
that are in some generic sense far away from equilibrium, such as driven diffusive systems [60], sheared 
fluids [61] or gelled [62] suspensions of Brownian particles, a living cell [63] or RNA molecules or protein 
domains unfolding upon application of an external pulling force [64], etc. Can, for such systems, still gen- 
eral statements be made that save us from having to solve the full dynamics? One approach tries to extend 
the classification scheme of universality classes [46] based on the renormalization group and equilibrium 
critical points to systems far from equilibrium [65] (see section 4). Closer related to the topics so far dealt 
with in the present section are two interesting lines of investigation that try to generalize the fluctuation 
dissipation theorem, which we recall is a generalization of the Stokes-Einstein and Green-Kubo relations. 
Both are to date still controversial and the subject of heated debates (e. g. [66]). 

The first is concerned with the intensely discussed paradigm of glassy dynamics [67]. Glassy dynamics 
is characteristic not only of structural glasses (such as colloidal, polymeric, molecular, or metallic glasses) 
and spin glasses but of a much wider class including many soft [68], disordered [69], granular [70], and in 
particular biological (e. g. proteins [71]) systems. As generic explanation for the origin of glassy dynamics 
several tentative mechanisms have been proposed (purely kinetic models, drop models, trap models, free- 
energy landscapes, etc. [68, 72-74]), and their further exploration to unveil the intricate roles played by dis- 
order and fluctuations remains a formidable challenge in theoretical physics. In contrast to the dynamics in 
supercooled liquids that approach the glass transition from the liquid side, which is rather well understood 
(see our discussion towards the end of section 6), the emphasis here is on the non-equilibrium dynamics 
after falling out of equilibrium, which is characterized by aging (i. e. relaxing indefinitely, in general in a 
non-universal manner, without ever reaching equilibrium) [73]. Yet, a generalized fluctuation-dissipation 
theorem (t > t') 



(J) 
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may still hold [73] (see figure 2). Close to equilibrium, T is the ordinary ambient temperature T and both 
sides depend on the time difference t — t' only. Equation (7) is then just a time-dependent generalization 
of Eqs. (2), (5) and the Stokes law F = (v. the ordinary equilibrium fluctuation-dissipation theorem. Far 
from equilibrium (e. g. below the glass transition temperature), this is still the case for t — t' <C t', where 
t' is identified with the waiting time after falling out of equilibrium. For t — t' ^> t' time-translation 
invariance and the relation between T and the ambient temperature are lost. The response and fluctuation 
functions in Eq. (7) are then dominated by their aging parts. The hope is that the effective temperature 
Tit') will still exist and be independent of the variable (here x) and consistent with alternative definitions 
of effective temperatures that are a priori independent of Eq. (7) 9 . This, in fact, would justify the notion of 
a generalized or quasi fluctuation-dissipation theorem [75]. While the theoretical understanding of these 
issues emerged from mean-field spin glass models [76], a corresponding development for structural glasses 
(e. g. suspensions of Brownian particles) is still at its beginning [77]. 

As a second line of attack towards the realm far from equilibrium, we want to mention the so-called 
non-equilibrium work relations ox fluctuation theorems. These are theorems supposed to remain valid 
beyond linear response, far from equilibrium, and in particular also for small systems 10 , away from the 
thermodynamic limit, where transient violations of the second law of thermodynamics become increasingly 
likely. The fluctuation theorems precisely say how likely. The recent excitement about these developments 
was initiated by the chaotic hypothesis of Gallavotti and Cohen [79] and independent work by Jarzynski 
[80], and it has even been suggested that non-equilibrium statistical mechanics is presently undergoing a 
foundational period comparable to that in equilibrium statistical mechanics a century ago [81]. The field 
has last been reviewed in Refs. [82-86]. Here, we try to sketch the basic idea. Consider a system prepared 
in thermal equilibrium with its environment at temperature T (e. g. an RNA molecule). The system shall 
be repeatedly subjected to an identical perturbation protocol such that the external work W is performed 
on it. In contrast to the reversible work Wq, which is a state function, W fluctuates from one repetition to 
the next. If the average over the (infinitely many) realizations of this process is denoted by an overbar and 
the dissipative work by W = W — Wo, Jarzynski's equality may be written in the form [85] 



e W/k B T = i ( 8 ) 

Obviously, the dissipative work must sometimes be positive as well as negative; i. e. sometimes mechanical 
work will be gained from the heat bath (see figure 1). This is indeed what Einstein's theory of fluctuations 
[38] requires. Already if the protocol does not take the system too far from equilibrium, so that fluctuations 
of W can be assumed to be Gaussian distributed about an average dissipation W, there will be some finite 
weight in the tail with negative W. This is not only consistent with the second law of thermodynamics, the 
latter is in fact immediately deduced from Eq. (8) via Jensen's inequality (e^ > e x ). Equation (8) can be 
used to infer the reversible work Wo from many repeated measurements of the irreversible work W. This 
has been tested experimentally, e. g. for a Brownian particle dragged through water [87, 88], and for RNA 
[64] and a number of other biomolecules subjected to mechanical forces (see Ref. [86] for an overview and 
suggestions for further experiments). 

4 Stochastic processes: the universal toolbox 

Universality: the reign of large numbers 

In section 2 we praised the elegance of Einstein's derivation compared to Smoluchowski's attempt to con- 
struct a microscopic dynamical theory of Brownian motion. This is not merely a matter of elegance, 



The effective temperature may, however, exhibit a more complicated time dependence of the form T({x(t)x(t'))), as observed 
in certain spin glasses, or simply diverge (corresponding to a universal "trivial" aging dynamics), as e. g. in ordinary coarsening 
during an equilibrium phase transition. 

Generally, care needs to be taken that the notion of temperature is still valid, which can be problematic for small systems [78]. 
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though. While the strength of Smoluchowski's approach lies in its ability to give insights into the mi- 
croscopic mechanism of Brownian motion, Einstein's approach relied on the idea that certain aspects 11 of 
fluctuations are independent of the details of the underlying microscopic dynamic processes, as long as 
they are sufficiently many and sufficiently uncorrected in time and space (both aspects also played a role 
in section 3). If this is so, can one develop a universally applicable calculus of stochastic phenomena? 
Indeed, this is essentially what the theory of stochastic processes aims at. 

Einstein's derivation utilized the Chapman-Kolmogorov equation Eq. (3) that translates the rules of a 
postulated simplistic time-discrete jump pr ocess of individual Brownian particles encoded in the jump 
propabilities tp T (Q, to the particle concentration n(x, t). Equivalently he could also have started from a 
formulation that is discretized in space and continuous in time and write down a master equation 

j t P{x, t) = Y}w^ x P(t t) - w x ^P(x, t)} (9) 

for a single particle with probabilities P(x, t) = n(x, i)/N and jump rates w^ x . Upon taking the limit 
of a large number of independent jumps, one obtains along either way a coarse-grained description in 
form of the diffusion equation Eq. (4). It is important to note that in taking this limit, some information 
about the underlying stochastic process, e. g. about the trajectory x(t), is irretrievably lost. The Chapman- 
Kolmogorov equation Eq. (3) and the master equation Eq. (9) provide a more detailed or microscopic 
description than the diffusion equation, Eq. (4). Indeed, Einstein exploits this fact by implicitely assuming 
that the difference between the actual microscopic dynamics, which Smoluchowski struggled to calculate, 
and the simplistic random walk will not matter for the end results Eqs. (4), (5). The idea becomes maybe 
most transparent from yet another equivalent approach Einstein could have taken. Starting from the discrete 
trajectory 

N 

z(i) = ^£(nr) ( 10 > 

71=1 

of a single particle with independent and identically distributed random jump lengths £, he could have 
appealed to the law of large numbers. Provided the first two moments of the jump distribution <p T (£) 
exist, the cental limit theorem guarantees that the limiting probability density for the individual particle 
trajectories x(t) starting at t = at the origin is given by the Gaussian [91] 



P(x, t) = , exp 
V^Dt 



x 2 
ADt 



(11) 



For the independent Brownian particles Einstein considered, this statement is of course equivalent to saying 
that the particle concentration P(x, t) obeys the diffusion equation, Eq. (4). In fact, Eq. (1 1) is nothing but 
the Green function of Eq. (4), i. e. its solution to the initial condition P(x, 0) = 6(x). 

Given such a powerful attractor of probability distributions, it is an intriguing question, whether there is 
a particularly simple member of the class of all stochastic processes with probability densities that converge 
to Eq. (11). Is there a simplest representative trajectory corresponding to the diffusion equation? This 
question was positively answered by Wiener in 1921 [92]. Inserting Eq. (11) with a short time interval 
t = t for the jump probabilities <p T (£) in the Chapman-Kolmogorov equation Eq. (3) and iterating the 
latter one arrives at the expression 



P{x,t) = [ jf-^U e-Ef.-o 1 ^^ _ ( X{t) V[x{r)\ e-^^mf^ ([2) 
J 7Li V^itDt Jx(o)=o 



^ \ZAitDt Jx(o)=o 



11 Precisely which aspects, is currently still an intriguing open question (see e. g. [89, 90]). 
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The final path-integral representation is obtained upon taking the continuum limit N — > oo with Nt = t, 
and has a well-defined meaning within the framework of measure theory [92]. Wiener could also char- 
acterize the typical paths corresponding to the limit taken in Eq. (12) and show that they are unphysical 
fractals, continuous but non-differentiable curves (see Fig. 4). One can imagine them to be constructed 
by a self-similar continuation of the large-scale structure of a trajectory of a Brownian particle to arbi- 
trarily short distances. The intricate details of the actual physical process at short scales together with 
the corresponding characteristic microscopic length scale are simply wiped off to give way to the ultimate 
simplicity of ideal self-similarity. The innocent looking prescription given after Eq. (12) actually encodes 
a renormalization of the original discrete stochastic model to a mathematically more convenient scale-free 
phenomenological model, known as the Wiener process. 



space 





Fig. 3 Illustration of the concept of a path in- 
tegral: trajectory of a Brownian particle starting 
at position xo at time to and passing through N 
gates of size dxi, dx2, ■ ■ ■ ,dxN-i,dx at times 

tl,t%, ■ • • , tjV-l, t. 



Fig. 4 Five different realizations of a random walk 
on a one-dimensional lattice illustrating the fractal 
nature of the trajectories, which in the continuum 
limit are continuous but non-differentiable. 



The art of mapping: from Brownian motion to quantum mechanics and back 

From what was just said, it is natural to expect that the Wiener process finds applications in many different 
areas far beyond the proper domain of Brownian motion, which is indeed the case [93]. Here we restrict 
ourselves to a few examples. Among physicists even better known than the Wiener measure Eq. (12) is 
the Feynman path integral [36], which expresses the propagator or Green function of a non-relativistic 
free particle in quantum mechanics as a path integral as in Eq. (12) but with an imaginary exponent. More 
precisely, the prefactor —1/2D is replaced by Within the path integral picture one can interprete 
the uncertainty relation of quantum mechanics as the analog of the diverging apparent velocity of the 
Brownian particles that worried experimentalists before 1905. Note that the exponent in Eq. (12) has the 
form of the Lagrangian of a classical free particle. This analogy can be extended [94] to include a potential 
energy —U(x) under the integral in the exponent of the Wiener measure. The corresponding Feynman 
path integral (the Feynman-Kac formula) can intuitively be understood as an analytic continuation of this 
extended Wiener measure in the complex plane, but a proper mathematical formulation is more subtle 
[93]. The diffusion (or Fokker-Planck) equation corresponding to the extended Wiener measure is Eq. (4) 
extended by a term —U(x)n(x,t), corresponding to an annihilation probability. The diffusion equation 
corresponding to the Feynman-Kac propagator is of course the Schrddinger equation (in imaginary time). 

Well known to physicists, this workhorse of quantum mechanics was the starting point for the develop- 
ment of various analogies. Building on the pioneering works, many important applications of path integrals 
have been developed in statistical physics and quantum and classical field theories, as well as in polymer 
physics and finance (where quantum mechanics was essentially mapped back onto Brownian motion), each 
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field contributing its own characteristic difficulties to the basic formulation. Out of the innocent jittering of 
a pollen grain has grown a universal and powerful tool of theoretical physics, which has helped to unravel 
intimate connections between what seem to be vastly different phenomena at first sight. 

We briefly sketch some key developments in the application to static conformations of flexible poly- 
mers. If segments of a flexible polymer could freely pass through each other, the Wiener process (in this 
context known as the Gaussian chain) would indeed provide an excellent coarse-grained description of the 
conformation if time t is identified with arclength s so that 6Dt is replaced by the mean-square end-to- 
end distance (R 2 ) of the polymer. In particular, the Stokes-Einstein relation is translated into the "rubber" 
force-extension relation 



d(R)i 
OF 



3k B T 



(13) 



which is easily recognized as the soft-matter analog of the Curie law for the inverse-temperature suscep- 
tibility of a paramagnet. Both become "stiffer" with increasing temperature, due to the disarranging action 
of thermal fluctuations. However, except under so-called 8— conditions, where the self-avoidance of the 
chain is screened by an effective self-attraction, the simple random walk is unrealistic because of this ide- 
alization. Real polymers conformations cannot intersect and are characterized by self-avoiding random 
walks. The Edwards model [95, 96] represents self-avoidance by substituting the potential in the extended 
Wiener measure by a singular repulsion U(r) oc 6[r(s) — r(s')]. This model has been the starting point 
of much theoretical work in polymer physics [97] and has been generalized to other self-avoiding man- 
ifolds [98]. It is amenable to a rigorous treatment with perturbation series and renormalization methods 
based on the path integral [97,99]. Notably, as has been pointed out by de Gennes [95], the problem of a 
self-avoiding random walk is related to critical phenomena, as it can also be mapped onto the n — > limit 
of a 0(n)— field theory. Figure 5 shows snapshots of fluorescently labeled DNA molecules adhering to a 
two-dimensional lipid membrane [100]. Measurements of the radius of gyration R g of the polymer have 
very well confirmed the prediction R g oc L 3 / 4 for the ("super-diffusive") conformation of a self-avoiding 
random walk in two dimensions. However, the interest of physicists for DNA is by no means limited to 
such very coarse and generic properties. Over the last 10 years DNA has become a favorite model system 
for a whole new field of investigation one might call "bio-polymer physics" [101-104]. 



> 1 «l > *? 



Fig. 5 Snapshots of fluorescently labeled DNA molecules adhering to a two-dimensional lipid membrane. Courtesy 
of J. Radler and B. Maier [100]. 



Self-similar renormalized models are mathematically convenient but fail to indicate where the simpli- 
fied universal representation breaks down. If one needs to incorporate more details, one has to resort to 
the level of the Chapman-Kolmogorov equation, Eq. (3), or master equation, Eq. (9). The latter has many 
features in common with the second quantized Schrodinger equation for many-particle quantum systems. 
In particular, both are linear differential equations of first order in time. Since each of the possible states 
a of the system can be labeled by a set of occupation numbers n(x, t) — > {rii(t)} (e. g. of lattice sites, 
energy states, chemical states, etc.), one is immediately led to a formulation in terms of an operator algebra 
in analogy with the second quantized formulation of quantum mechanics (for a short review and a list of 
references see e. g. [105]). Upon introducing boson ladder operators with the standard commutation rela- 
tion [en, Oj] = Sij, one can define Fock states |a) = Ili( a i |0) y i a act ing with the ladder operators on 
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the "vacuum" state |0). Then the master equation for the probability P({m}, t) of finding the system in 
the state with {n^} particles on site {i}, respectively, is equivalent to a Schrodinger equation in imaginary 
time 

^|*) = (14) 

for the state vector \^f) = ^2 a P({rii}, t)\a). Of course, there are also differences from quantum mechan- 
ics. Most importantly, the "Hamiltonian" TL is in general non-Hermitian. Still, most if not all many body 
techniques can be applied. In particluar, it is possible to transform the second quantized Hamiltonian into 
a field theory upon using coherent-state path integrals [65]. 

The similarities of the independently developed formalisms have led to a successful description of 
classical statistical systems far from equilibrium [105]. They not only allow for an application of field 
theoretical methods but also for a fruitful exchange of results and concepts between seemingly unrelated 
phenomena in quantum and classical physics. Very recently, it was e. g. suggested to represent gene 
expression as a many-body quantum problem [106], which leads to such unexpected analogies as that 
between the stochastic dynamics of a single gene switch and the quantum mechanics of the spin-boson 
model of a two-site polaron. 



5 A noisy world: the ubiquity of stochastic processes in biology 

Genetic drift and the theory of neutral evolution 

Today we live in a "noisy world" where it is difficult to find a quiet place to contemplate. Noise in this 
context is simply unwanted sound and certainly deleterious. However, at the very origin of our lives are 
random processes, i. e. noise. Luckily, one should say, because otherwise we would largely be deprived of 
our individuality and doomed to live in a world of clones. This insight we owe to Mendel, who deduced 
from his experiments with pea flowers the existence of genes present as pairs of alternative forms (called 
alleles) in each (somatic) cell. He assumed that reproduction happens in a two step process, where equal 
separation of these gene pairs into the gametes (egg and sperm) is followed by random fertilization, i. e. 
random union of one gamete from each parent to form the first cell (zygote) of a new progeny's individual. 
This process is random precisely because there is Brownian motion on a molecular scale. 

Starting at about the same time as Einstein's fundamental paper on Brownian motion the re-discovery 
of the Mendelian laws prompted the emergence of a new discipline, now known as population genetics 
(for an historic account see e. g. Ref. [107]). Pretty early on it was realized that random fertilization 
in finite populations leads to a phenomenon called "genetic drift" [108]. To illustrate this point let us 
look at the following toy model of population genetics (see Fig. 6), where one considers only one gene 
locus in a population of diploid 12 organisms that has two alleles A and a. We assume that each individual 
releases the same number of gametes, so that the allelic frequency in the population of gametes and the 
parental population are identical. From this large pool of gametes the next progeny is obtained by randomly 
chosing N pairs of gametes ("random mating" or "panmixia"). In the language of sochastic processes this 
is a sequence of 2N independent Bernoulli trials. If we denote by x the relative frequency (probability) 
for the allele A in the parental population, the probability distribution for the number n of alleles A in the 
next generation is given by the binomial distribution 

B{n,N)= Nl - x n (l-x) N - n . (15) 
n\(N — n)\ 

Therefore, even if there are no deterministic evolutionary forces acting on a population (such as natural 
selection), sampling errors during the production of zygotes may cause a change in the frequency on a 

12 The genome is made up of one or more extremely long molecules of DNA that are organized into chromosomes. The body 
cells of most plants and animals contain two genomes. These organisms are called diploid. For more details the reader is 
referred to a recent edition of some textbook in genetics, e. g. [109]. 
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particular allele in a population; this phenomenon is called genetic drift. As can be read off from the 
variance of the binomial distribution, D{x) = x(l — x)/2N, genetic drift is expected to be large in small 
populations and vanishes for infinite population size. In the language of Brownian motion "genetic drift" 
is actually no drift at all, but a random walk of the stochastic variable x (gene frequency). Two examples 
of such walks either leading to extinction or fixation are shown in Fig. 7. 



first generation 




random sampling 



male gametes female gametes 
o • o o 

o o • • 



fertilization 




second generation 



Fig. 6 Panmixia: a toy model of population genetics. 
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Fig. 7 Genetic drift. Some representative sample 
paths for random gene frequency drift. 



In complete analogy to Einstein's formulation of Brownian motion one can now consider an ensemble of 
populations and derive a differential equation for the conditional probability density n(x,p; t) that the gene 
frequencs lies between x and x + dx at time t, given that the initial gene frequency was p. It is interesting 
to note that such a formulation was given by Fisher [108] without being aware of parallel developments in 
physics. The corresponding Fokker-Planck equation reads 



with D(x) = x(l - x)/2N, 



(16) 



which can be interpreted as a diffusion process with progressively slower motion near the boundaries of 
gene frequency space, i. e. x = and x = 1. Similar models with spatially varying diffusion coefficients 
have more recently been studied in the context of the voter model [1 10, 1 1 1] and of some catalysis reaction 
models at surfaces [112]. 

An essential feature of the random "drift" in gene frequency space is that there are boundaries. Once 
a given allele A reaches x = it becomes extinct and the random process stops. Similarily if the allele 
A reaches a frequency of x = 1 it is said to be fixed (or equivalently allele a becomes extinct). In other 
words, what we are facing here is a diffusion process with two absorbing boundaries. One may now ask 
the following questions. What is the time evolution for the distribution of unfixed classes with < x < 1? 
How likely is it for an allele to become fixed and what is the statistics of the number of generations until 
extinction or fixation of a certain allele occurs? Both of them can be answered by analyzing the Fokker- 
Planck equation Eq. 16, whose application to population genetics has been pioneered by Kimura, one of 
the founding fathers of the neutral theory of evolution 13 . He showed [113] that in the asympotic limit 
of a large number of generations the distribution function of the unfixed classes becomes flat and decays 
exponentially as n(x,p;t) — 6x(l — x) exp(— t/2N). For answering the second question there is an 
alternative route using concepts from game theory. To see this let us shortly digress to a phenomenon 
called "gambler's ruin". Say the game simply consists of repeatedly throwing a fair coin. Your initial 
capital is p (with < p < 1) and that of your opponent is 1 — p. The game is over when either one of 



The neutral theory of evolution asserts that the great majority of evolutionary changes at the molecular level are caused not by 
Darwinian selecion but by random drift of selectively neutral mutations [107]. 
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you goes broke. What are your chances to win the game? The answer is simple (see e. g. Ref. [114]). 
Your chances are p, i. e. proportional to your initial capital, which can be rather low if you try to play 
against someone with large reserves 14 . This is exactly the game played by the genes in the above toy 
model of population genetics. Hence the probability that the allele A eventually becomes fixed ("wins") is 
identical to its initial frequency p. But how long will it take to win? This is a first passage time problem 
(for a review see e. g. Ref. [115]) and can be solved by using a cousin of the Fokker-Planck equation, the 
backward Kolmogorov equation [116]. If one starts with just a single "invading" allele A with a frequency 
of p = 1/2N in a population of size N, one finds [117] that both the mean and the standard deviation of 
the time to fixation grow linearly with the size of the population. 

This simple toy model of population genetics with one gene locus, two alleles and random mating is, of 
course, much to naive to account for the variety and complexity of life forms on our planet. On a molecular 
level, there are complex changes in DNA sequence, which may either be random-like point mutations and 
frameshifts or more systematic cut-and-paste mechanisms like insertion sequences and transposons [118]. 
On the phenotype level, there are interactions among the organisms in a population and with the sur- 
rounding environment, which may itself be dynamic and noisy. The fitness, i. e. the number of expected 
offsprings, of an individual is a (highly non-linear) function of such interactions. The way from genotype 
to a phenotype is paved by networks of coupled biochemical reactions and feedback signals. The list goes 
on and would lead us considerable off track; for more information the interested reader is recommended to 
consult some recent introductory review articles [119, 120]. 



Stochasticity in chemical reactions: the reign of small numbers 

The machinery of biological cells consists of networks of molecules interacting with each other in a com- 
plex manner. Many of these interactions can be described as chemical reactions, where the intricate pro- 
cesses which occur during the encounter of two molecules are reduced to reaction rates 15 . For example, 
an autocatalytic reaction for an enzyme E is modeled as a probabilistic self -replication process E — » 2E, 
which occurs at a constant probability per unit time (reaction rate) k. For high concentrations of the 
enzyme, a standard approach in chemistry is to write down rate equations for the time evolution of the 
average concentration n(t) 

Such an approach evidently assumes that the time evolution of a chemically reacting system is both deter- 
ministic and continuous. However, the actual process is neither of both since molecular population levels 
change stochastically and by discrete integer amounts. 

Obviously, a deterministic description is applicable only if deviations from the average number of 
molecules are negligibly small. As noted long ago by Delbriick [123] the validity of this assumption 
breaks down in small chemical systems, where the concentrations of reacting species are extremly low and 
the statistical fluctuations of atomism can not be avoided. Then, instead of the average concentration n 
one now has to consider the full probability distribution function P n (t) of the number of enzymes n at 
a given time t. This function plays the same role as the particle density n{x, t) in Einstein's treatment 
of Brownian motion or the probability density of gene frequency n(x,p, t) in our discussion of genetic 
drift. Instead of continuous random variables (particle position, gene frequency) we are now dealing with 
the discrete random variable n(t), the number of enzymes at time t. The time evolution for P n (t) can 
immediately be written down if one pictures an ensemble of enzymatic reactions, where P n is the rela- 
tive frequency of samples in the ensemble with exactly n molecules. The frequency P n (population size) 

14 But history teaches us that even the big players with large reserves may sometimes stumble and fall [7]. 

^ Reaction rate theory was pioneering by Kramers [121]; for a recent review see e. g. [122]. He considered a simplified model 

system where a Brownian particle is trapped in a one-dimensional well representing the reactant state which is separated by a 

barrier of finite height from a deeper well signifying the product state. 
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of this particular class is reduced by each self -replication process within the class and enhanced by each 
self-replication processes in the neighboring class. Then the time evolution of P n {t) can be written as a 
"gain-loss" equation 

j t P n (t) = (n - l)fcP n -i(t) - nkP n (t) , (18) 

a specific example of the master equation introduced in section 4. 

How are the results of a stochastic description, Eq. 18, different from a deterministic approach, Eq. 17? 
The most striking one is, of course, that despite the exponential growth, n(t) = ra exp(fct), promised by 
a deterministic description there is a chance to loose it all, n{T) — (gambler's ruin). Other features of 
the stochasticity become evident upon solving the master equation for reactions initiated by one particle. 
Then the probability distribution for large times t has the asymptotic form 

P r ](t) = 1 e~ n ' n with n = e kt , (19) 

which is by no means normal, but exponential. Note that at all times the probability of finding a single 
particle is greater than the probability of finding any other specified number of particles. In addition, the 
standard deviation is idential to the mean of the distribution. Delbriick also points out that there are strong 
fluctuations in the waiting time for attaining a certain threshold amount of reaction products [123]. 

An even more dramatic illustration for the importance of stochasticity in chemical reactions is given by 
the pair annihilation process A + A — ► in one-dimensional systems. The rate equation, dn/dt = —kn 2 
would predict an algebraic decay as n(t) <~ i -1 . However, the actual asymptotic decay is much slower, 
n(t) ~ i -1 / 2 . This slow decay can be traced back to the reentrance property of random walks in less 
than two dimensions, which implies that within a volume of order \f~Di, all particles are annihilated with 
probability one [124]. A more sophisticeted line of argument uses methods adapted from quantum field 
theory [125]. The pair annihilation process is the simplest member of a broad class of models known as 
diffusion-reaction models exhibiting absorbing states, recently reviewed in Ref. [126]. That such systems 
show unusual dynamics becomes obvious by asking the innocent looking question: How does a system 
evolve towards its steady state? Sometimes the answer is quite simple and the relaxation process is merely 
an exponential decay. If the deviations from equilibrium are small Onsager's regression hypothesis [127] 
asserts that the relaxation is governed by the same laws as the fluctuations in equilibrium. But this hypoth- 
esis certainly fails for systems with an absorbing steady state. Here there are no fluctuations in the steady 
state but the approach towards the absorbing state is critical in the sense that it exhibits slow power-law 
decay and universal scaling behavior [128]. 

Before getting too much sidetracked let us come back to biology. Here it has been emphasized recently 
that living cells possess very low copy numbers of many components, including DNA and important regu- 
latory molecules [129]. As a consequence, stochastic effects such as those illustrated above may play a key 
role in biological processes. A quite striking example are observations [130, 131] that clonal populations 
of cells (i. e. genetically identical cells) exhibit substantial phenotypic variation. Recent experiments even 
allowed to distinguish between the stochastic noise inherent in the chemical reactions controlling gene 
regulation and fluctuations in the number of other cellular components [132]. These observations are only 
the tip of an iceberg. "Noisy" mechanisms may well be rather the rule than the exception in the molecular 
world of a cell [133]. 

6 The Brownian particle: probe and effigy of the intangible 

Probing molecular fluctuations with Brownian particles: "the buoy in the surge" 

Precisely in the way proposed by Einstein [2], Brownian particles have ever since been used as probes for 
the molecular motion of the surrounding medium: recently in particular for probing the local viscoelastic 
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response of soft, often biological, materials, where elastic moduli are of the order of fc^T or even due to 
ksT (i. e. of entropic origin). A zoo of methods based on Brownian motion of tracer particles has been 
invented: (z) particle-tracking techniques 16 that record the path of one particle [17, 135, 136] or the corre- 
lated fluctuations of two particles [137, 138]; (ii) diffusing light spectroscopy [139, 140], which iterates our 
basic theme by exploiting not only the diffusion of colloids but also the diffusion of light they scatter; (Hi) 
Fourier-imaging- [141] and fluorescence- [142] correlation-spectroscopy, which are more similar to clas- 
sical dynamic scattering techniques; etc. According to the fundamental connection between fluctuations 
and response provided by Eqs. (2), (5), (6) the fluctuations measured by these passive techniques contain 
the same information as obtained by active microrheological linear-response measurements (e. g. [143]). 
One can roughly classify the passive methods into those working in Fourier space ("scattering") and those 
working in real space ("microscopy"). 

The basic principle of scattering methods probing Brownian dynamics is well illustrated by the example 
of particle sizing, which is a widely applied standard method. One shines a laser beam of wavelength A at 
a dilute suspension of N identical particles in a viscous medium and records the light intensity scattered 
at an angle 9 to the incident beam. The auto-correlations of the scattered light then decay proportional to 
(the square of) the dynamic structure factor 

S q (t) = — ^( e «J-M*)-' r 3 ( )]) di ^ t0 — ^( e ^-^.(t)) = e -i 2 (Sr t (t) 2 )/6 = e -q 2 Dt _ ^0) 

ij i 

Here, q is the scattering vector, q = \q\ = Ann sin(0/2)/A with n the refractive index, and we used that 
dilute particles are uncorrected. Applying Einstein's Eq. (5) together with the Stokes relation ( = Qn-qa 
for the friction of a spherical particles in a viscous solvent, one can determine the diameter a of the particle 
or alternatively the vicosity r\ of the solvent. 



Brownian particles in complex fluids: "the buoy in the pudding" 

Although the idea of exploiting such kind of methodology to measure material properties in small sam- 
ples and even inside living cells is also a rather old one [144], its applications mulitplied over the last two 
decades or so with surgeing interest in highly disordered, soft, and/or biological matter. Reports of sub- 
stantial deviations from Eq. (5) in such systems abound. A common observation is that of superdiffusive 
(a > 1) or subdiffusive (a < 1) dynamics, which was studied theoretically in great detail [72], and which 
is characterized by 

(Sr 2 )(xt a witha^l. (21) 

Such "strange kinetics" [145] has even prompted the development of a dedicated "fractional calculus" 
[146-148]. It appears in many areas of science (nonlinear dynamics, growth phenomena, fluctuating man- 
ifolds, etc.). A simple way to generate such anomalous diffusion is by assuming jump rates </v(0 in 
Eqs. (3) (10) that decay too slowly with £ (i. e. like 1 < p, < 3) for the ordinary central limit theo- 
rem to apply, or, equivalently, for constant step length by a broad waiting time distribution. For Brownian 
particles, it naturally arises as intermediate asymptotics whenever the solvent is not a simple liquid but a 
complex (viscoelastic) fluid with some internal structure that exhibits long-range spatial or temporal cor- 
relations; i. e. it has more to do with a complicated material property of the solvent than with Brownian 
motion as such, and thus leads us slightly off-topic. Biological cells certainly provide us with ultimate 
complexity, which makes deviations from Eqs. (5), (20) anything but surprising. It is therefore very sug- 
gestive to interprete subdiffusion in cells as a consequence of cytoplasmic heterogeneities and crowding 
[149, 150], but attempts to explore the actual connection are still limited to model systems [151, 152]. If 

16 Roughly speaking, even the latest computer-aided confocal video-microscopy setup for particle tracking is but a costly late 
descendant of the technique pioneered by Henri in 1908 [134], when he first put a film camera onto a microscope to record the 
wild paths of Brownian particles in order to scrutinize Einstein's prediction, Eq. (5). 
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not due to complex (biological) organization, one generally expects long-range correlations to be self- 
similar, arising e. g. as a consequence of a closeby equilibrium critical point, turbulence, or some other 
non-equilibrium structure. A nice example of temporal correlations is provided by simple shear flow: it 
gives rise to a peculiar cooperation of statistical mechanics and fluid mechanics known as Taylor diffu- 
sion [29, 153], which effectively speeds up the diffusion of (dilute) suspended Brownian particles. Stirring 
your coffee thus helps twice to distribute the suger, via turbulent convection on large scales and via the 
induced superdiffusion on shorter scales. A common paradigm based on spatial correlations is diffusion 
in a porous medium with a very broad pore size distribution, idealized as diffusion on a fractal [154, 155]. 
More intriguing from a theoretical perspective is the case, where subdiffusion arises without any non- 
trivial structures in configuration space as a consequence of a complicated phase-space structure. This 
kind of "dynamically self-generated" subdiffusion is often characterized as glassy dynamics (see section 3 
and below). 

Long-time tails: "the buoy's bow wave" 

Maybe most suprisingly, already the dynamics of a single Brownian particle in a simple Newtonian solvent 
has some subtle properties not anticipated by Einstein and Perrin. In his derivation, sketched above, Ein- 
stein implicitely assumed <p T (Q to have properties such that the expansion of the Chapman-Kolmogorov 
equation, Eq. (3), nowadays known as Kramers-Moyal expansion, can be truncated after the V 2 — term, 
thus yielding a simple hydrodynamic equation, the diffusion equation, Eq. (4). Such conditions are math- 
ematically realizable [27,28, 156] but not necessarily physically relevant [116, 157]. They amount to a 
description of Brownian motion as a Gaussian Markov process (a notion that did not yet exist in 1905). 
Although this may seem to be a very plausible assumption given the law of large numbers and the strong 
scale separation between the Brownian particle and the molecular scale, it only holds to a first approxi- 
mation, or for t — > oo. Hydrodynamic conservation laws, as e. g. contained in the Stokes equation for 
the solvent, induce long-range auto-correlations that effectively retard the approach to the Markov limit. 
At any finite time, they contribute transient corrections, so-called long-time tails, to Eq. (4) and more- 
over render Eq. (3) and the subsequent expansion ill defined [158-160] (for overviews see e. g. [157, 161] 
and references therein). A quick derivation exploits Eq. (6). Recalling that the dynamics of the solvent 
is itself governed by a diffusion equation (see footnote 2, section 2, above) we conclude that the mo- 
mentum of the Brownian particle is diffusively dissipated, hence shared with a growing solvent volume 
A(t) 3 oc (^t) 3 / 2 , so that the velocity auto-corelations in Eq. (6) decay like t~ 3 / 2 . The effective diffusion 
coefficient (Sr 2 (t))/6t at finite time t is obtained by setting the upper limit of integration in Eq. (6) to t, 
which immediately gives [157, 162] 



where 9 is a hydrodynamic time scale that depends on £ and the particle/solvent mass ratio. Thus, normal 
Brownian motion actually is slightly faster than diffusive, Einstein's formula is only slowly approached, 
and there is no simple general way to ammend the diffusion equation to take this into account. The some- 
what uncomfortable picture [157], which first emerged from calculations for dense gases based on the 
Boltzmann equation [158], has sometimes been summarized as Dorfman's lemma [163]: all relevant fluxes 
are nonanalytic functions of all relevant variables. Neither does a gradient expansion of the diffusion cur- 
rent exist (the Burnett-coefficients appearing in the expansionn are found to be infinite), nor do density 
expansions exist of the diffusion coefficient and other transport coefficients as one might expect in analogy 
with the virial expansion of thermodynamics. Luckily, the corrections in Eq. (22), though of fundamental 
interest, are in general numerically too small to spoil seriously any of the numerous applications mentioned 
above. We note in passing that related problems arise in statistical mechanics derivations a la Einstein and 
Smoluchowski of hydrodynamic boundary conditions (e. g. "no-influx" for a tagged Brownian particle in 
a colloidal suspension [164] or "no-slip" in a Newtonian fluid [165-167]). Again, typical microrheolog- 
ical applications may well get by without such subtleness [168]. However, one should be aware that tails 




(22) 
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do seriously affect Brownian motion in lower dimensions, leading to "hydrodynamic suicide" [51] in two 
dimensions. Even without considering solvent hydrodynamics, due to the otherwise subdominant long- 
range concentration patterns ("colloidal tails") implied by Eq. (4), Brownian motion in 1-dimensional 
pores exhibits anomalous, subdiffusive dynamics, so-called single-file diffusion [169]. 

Interacting Brownian particles: the upscaled microcosmos 

Apart from providing us with a large assortment of colloidal particles to probe the molecular world, the 
methods pioneered by Perrin, Svedberg, and others moreover opened the way to designing an effigy of the 
molecular world on the colloidal scale [170-172]. The natural laboratory of Gouy was successfully trans- 
formed into an artificial microcosmos [25], which obviously leads us slightly off-topic, again. Svedberg's 
ultracentrifuge helped to solve one of the big experimental challenges in colloid preparation: to prepare 
samples that have sufficiently well-defined properties, in particular a monodisperse size distribution to 
mimic as closely as possible the sameness of atoms. This is hard work: starting with one kilogram of 
gamboge, which is a gum resin named after Cambodia (from where it was originally imported), Perrin 
"obtained after several months of daily operations a fraction containing several decigrams of grains with a 
diameter of approximately three-quarters of a thousandth of a millimetre" [5]. But monodisperse micron- 
sized particles were a crucial prerequisite for the progress of colloid science. With an enormous toolbox at 
hand, chemists are now able to design almost any colloidal interactions one may wish for, far beyond what 
can be found among naturally occuring molecules. 

This gave experimentalists a formidable playground for testing highly idealized theoretical pet-models 
for many-body problems, such as the hard-sphere system and its various extensions. As anticipated by 
van der Waals in the 19th century, much of the packing structure of liquids can be understood from the 
hard-shpere like mutual repulsions between their constituent molecules [173]. Hard spheres also provide 
a beautifully simple example (paradigmatic for more complicated soft condensed matter systems [15] in- 
cluding living cells [174]) how macroscopic order emerges from microscopic chaos, namely the stochastic 
Brownian motion of the individual spheres [175]. Numerical simulations [176, 177] and experiments [178] 
have established that above a sharply defined volume fraction <j> — 0.494, Brownian fluctuations drive 
the particles into a crystal, because they gain space to wobble around when they arrange themselves on 
a (virtual) lattice 1 . Microgravity experiments [179] suggest that this crystal is made out of randomly 
stacked hexagonally close-packed planes. It becomes space filling at <f> = 0.545 and seems to develop a 
slight preference for face-centered cubic organization if it is compressed towards closest packing [180] at 
4> = 0.74. Colloidal crystallization has gained additional interest from those trying to crystallize proteins 
for structure analysis [181]. While, in this context, colloids are currently still mostly playing their usual 
role as visible model systems for globular proteins [182], programmable protein expression might in the 
future well overcome chemistry as the more powerful machinery for providing colloid physics with the 
ultimate design precision. 

If quickly quenched to volume fractions cj> > 0.58, hard spheres will provide yet another interesting and 
suprising feature closely related to our central theme: by their Brownian motion they trap themselves in 
an arrested non-equilibrium amorphous state that looks like a liquid but feels like a solid, and which may 
substantially delay or even impede crystallization: the colloidal glass 1& . From mode-coupling theory [185] 
a very characteristic (self-generated) subdiffusive behavior [186] during the critical slowing-down of the 
Brownian motion near this transition was predicted. This was beautifully confirmed (see e. g. Refs. [187, 
188]) for collective density fluctuations by light scattering [189], and for tagged particles by numerical 
simulations [190]. At least for hard-sphere like systems the theory can be regarded as legitimate — though, 



This illustrates the paradox that global order may help to increase individual freedom. 

To be fair to the photon [183], we note aside that not only diffusing Brownian particles may come to a halt, but also diffusing 
light. Recently, colloidal crystallization has gained some attention as a promising method to produce photonic crytals [184], 
which may become an attractive technology in the future (currently all the rage: hiding light under a bushel). 



adp header will be provided by the publisher 



19 



unfortunately, not completely rigorous — sucessor of Eq. (4) towards a more microscopic picture. It 
starts from a microscopic expression for the normalized dynamic structure factor <fr q (t) = S q (t) / S q (0), cf. 
Eq. (20), in terms of generalized fluctuating particle currents J 



For short times and q —* 0, this reduces to the Fourier transform of the diffusion equation. The kinetic 
coefficient D(q) generalizes the gradient diffusion coefficient D = D(q — > 0) to finite wave vectors. For 
long times, it becomes effectively reduced due to the term reminiscent of Eq. (6), through which the system 
remembers its past evolution (e. g. its long-time tails). To make practical use of Eq. (23), the microscopic 
memory kernel containing the Brownian fluctuations J has to be expressed in terms of the hydrodynamic 
variables <f>, which involves the uncontrolled mode-coupling approximation [187,192,193]. Yet, after 
iterative solution, Eq. (23) predicts in impressive detail how, upon increasing either the attraction or the 
volume fraction, the simple Brownian dynamics of a fluid suspension described by Eq. (4) gives way to 
slow glassy dynamics, and eventually to two disticint non-ergodicity transitions into two fundamentally 
different glassy states [194, 195]. 

7 Fluctuating manifolds: Brownian particles with soft internal modes 

Just as it is commonplace that history repeats itself, re-inventing the wheel is an everyday experience in sci- 
ence. So it could have been expected that researchers studying red blood cells under the miscrosope would 
be tempted to ascribe their flickering shape undulations to animate causes, thereby repeating the initial 
misinterpretation of Brownian motion as a sign of life by (most of) its early discoverers (see Ref. [196] for 
a discussion of the history). The story of Brownian motion was about to repeat itself for Brownian particles 
with soft internal degrees of freedom. Today, cell flickering is understood to be a physical phenomenon 
[196] that has an important biological implication, tough: the reduction of unfavorable adhesion via the 
so-called Helfrich repulsion [197] caused by the thermal undulations of the cell membrane. Similarly, the 
thermal conformational wiggling of the polymeric constituents of the cytoskeleton has been pinpointed as 
the origin of some characeristic viscoelastic behavior [168, 198, 199]. 

Brownian undulations . . . 

For fluctuating manifolds, not only the center of mass position or overall orientation but an (ideally) infi- 
nite number of internal modes are excited into a persistent dance by thermal forces. Brownian fluctuations 
tend to "crumple" low dimensional manifolds [98]. Their free and driven 20 conformational dynamics is the 
subject of a field of investigation that one might summarize by the term stochastic elastohydrodynamics. 
It is largely motivated by the current interest in soft and biological materials. Examples of fluctuating 
manifolds are flexible linear and branched polymers, microemulsions, foams, liquid crystals and numer- 
ous other soft surfaces, interfaces and membranes [96,98,204-208], but also vortex matter, line liquids, 
magnetic domain walls and the like [209, 210]. 

For the most natural case of an elastically uniform manifold of n = 1 or n = 2 dimensions [211], 
the force needed to excite internal modes will grow as a power law with the mode number q, equilibrium 
amplitudes a q and relaxation times will decrease as (a q ) oc q~^ and r q oc q~~*, respectively. For 7 > /?— n, 



The worst theory apart from all others that have been tried from time to time [191]. 

Opening the Pandora box of non-equilibrium one enters the vast topic of pinned or driven fluctuating manifolds [200-203] 
that comprises a wealth of complex non-equilibrium phenomena in many areas of condensed matter physics, which we will 
not pursue any further, here. 
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this is quickly translated into a subdiffusive' 



mean-square displacement 



(24) 



by transforming the sum into an integral. This is the appropriate generalization of Eq. (5) to real-space 
measurements of a small labelled patch on the manifold, which have indeed been very successfully applied 
to various biomolecules [212-214]. Similarly as for simple Brownian particles, Eq. (24) determines the 
decay of the dynamic structure factor for long times (t — ► oo), when only the self-correlations of patches 
of dimension q^ 1 matter, via Eq. (20). In contrast, for times shorter than the relaxation time of modes of 
wavelength q^ 1 , the fluctuating interference of light scattered from neighboring elements of the manifold 
within such patches has not yet averaged out and dominates the dynamic decorrelation of the scattered 
intensity. While the internal modes thus affect the decay law for long times, S q (t — ► oo) oc e~ Vq * , they 
merely renormalize the effective diffusion coefficient for short times, S q (t — > 0) oc e~ Vqt . Both limits 
encode information about the material properties of the manifold and the solvent in different combinations. 
(Which regime dominates the dynamic structure factor is determined by the elasticity of the manifold. For 
flexible manifolds the simple exponential behavior dominates, and for stiff ones the stretched exponential 
tail.) Under favorable conditions this can be exploited to probe these parameters [215] in analogy to the 
simple example of particle-sizing discussed above. Fluctuating manifolds thus add a new facet to the story 
of Brownian motion: subdiffusive behavior may not only arise due to complicated surroundings but equally 
well from (quite simple!) internal elastic degrees of freedom of the tracer. In particular in presence of exter- 
nal forces this opens up the possibility of rich non-equilibium Brownian dynamics even for such seemingly 
trivial objects as a thread in a Newtonian solvent [216]. If the long-range hydrodynamic self-interactions, 
mentioned above in connection with the long time tails, need to be included, they lead to considerable 
complications, as e. g. for flexible polymers [96] and for membranes [217], even in equilibrium. 

... in complex fluids 

Increasing the complexity by one more step, one can combine both extensions of simple Brownian motion, 
internal elastic modes and a complex solvent. Consider a long (tagged) polymer embedded into a dense 
melt or solution of other (identical) polymers. The question how the polymer diffuses through the medium 
is closely (but by no means trivially) related to the question how the whole medium is able to flow. Both 
are addressed by reptation theory [96,204,218], the prevailing phenomenological idea how this involved 
many-body problem can be disentangled and reduced to a description based on a single polymer laterally 
confined to an effective tube. The concept is popular for its intuitive elegance and its ability to rationalize 
data for the elastic modulus and the viscosity. One imagines the polymer as a random walk caught in a cage 
shaped like a curly tube of mean-square end-to-end distance (R 2 ) proportional to its length L. Taking 
into account that it has to diffuse back and forth all the way along the tube to completely disengage from 
it, and that the diffusion coefficient D\\ for diffusion along its own contour is inversely proportional to L, 
one deduces the disengagement time 



For shorter times, due to multiple inter-polymer collisions [99], semidilute solutions or melts respond 
elastically to shear, like a rubber with an elastic modulus of ksT per statistically independent segment of 
the tube. These segments (of length (R 2 )/L) act as its effective unit elements. At times longer than r r , 
entanglements can relax, and the whole medium can flow like a liquid. The large time r r therefore also fixes 
the ratio of the medium viscosity at long times to the shear modulus at shorter times, which implies that 
polymer melts and solutions are soft but tenacious, or "gloop". Support for the tube, intuitive starting point 



T r oc L 2 /D\\ oc L 3 . 



(25) 



If this inequality is violated, one does not obtain superdiffusion but the dynamics is dominated by the lowest mode. 
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of the reptation model, comes from direct visual observation of the Brownian motion of large biopolymers 
[219,220]. Towards a deeper understanding of how viscoelasticity emerges from the underlying basic 
laws of Brownian motion, microscopic approaches (theory [221] and simulations [222,223]) are currently 
pursued. 



Fig. 8 Direct visualization of the tubes around actin filaments in F-actin solutions. Courtesy of J. Kas [219]. 

Over the last decades, membranes and polymers have become the major companions of the simple 
Brownian particle as paradigms of soft condensed matter, which is what we ourselves and many of our 
natural and technological aides of everyday life (foods, cosmetics, detergents, plastics . . . ) are made of. 
Rather than continuing this list ad infinitum, we close this section by pointing out an ironic twist in the re- 
cent history of Brownian motion: with the physics of fluctuating polymers and membranes reasonably well 
understood, physicists have at last turned their attention back to the possible animate sources of Brownian 
fluctuations. The idea is to study the additional effect of "active noise", such as the undulations induced 
by molecular pumps, molecular motors, and switching channels embedded into a membrane [224,225], 
a theme that can be iterated in many directions; e. g. towards so-called active gels [226-230], solutions 
of biopolymers mixed with molecular motors, in which rich structure formation (asters and spirals, etc.) 
has been observed in vivo and in vitro. Interestingly, the "noise" produced by the biological activity was 
predicted to lead to a thousandfold increase of the amplitudes of the long-time tails [231]. 



8 Rectified Brownian motion: from "life force" to "living force" and back 

It is an amusing semantic detail from the history of science that kinetic energy was originally termed living 
force ("lebendige Kraft") by Helmholtz and his contemporaries, adopting Aristoteles' "vis viva". When 
Brown first observed the irregular motion of pollen granules, he was convinced that he had found animate 
objects driven by some kind of "life force" outside the realm of physics and chemistry, a vitalistic view 
still common at that time. He derived from Leeuwenhoek's terminology of "animalcules" (small animals) 
the name "molecules" for the rambling particles under his microscope. After their animal aura had been 
stripped off by the subsequent investigations, the term could be reused for those particles that carried the 
thermal living force then recognized as the mundane physical explanation of the hitherto mysterious life 
force. 

This section will deal with the origin of motion in living systems and will bring us back once again 
to the "living force" as the "life force" responsible for directed motion on a cellular scale. The origin of 
motion of living organisms was already debated in the 3rd century B.C., when the anatomist and physician 
Erasistratos of Ceos associated muscle motion with the "spiritus animalis". He imagined it as some kind 
of fluid or gas, which he termed "pneuma", that flows through hollow nerves as pipelines and makes the 
muscles swell and shorten. This pneumatic idea survived in the scientific community for a very long time. It 
was only with the invention of the microscope by van Leeuwenhoek (1632-1 723) that S wammerdam (1637- 
1680) was able to show that muscles contract at constant volume which invalidated pneumatic theories. 
Helmholtz in his famous paper about the law of energy conservation [232] may have been the first to 
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emphasize that the mechanical energy produced by living organisms is just transduced chemical energy . 
These ideas paved the way for a physical understanding of muscle contraction. A discussion of the early 
work on theoretical models, which already contains some of the ideas of more recent work discussed below, 
can be found in reviews of A.F. Huxley [233] and Hill [234]. 

Modern experimental techniques [235] have lifted the veils and allowed us to look at the causes of 
biological motion on a molecular scale. This resulted in the following picture [236]. Biological mo- 
tion is caused by (nanometer sized) motor proteins, a highly specialized class of enzymes which are able 
to transduce the energy excess in the chemical hydrolysis reaction of ATP (adenosine-triphosphate) into 
mechanical work. An important subclass are cytoskeletal motors, which by sequence similarity are clas- 
sified into three families (myosins, kinesins and dyneins). These motors "walk" along one-dimensional 
molecular tracks consisting of quite stiff protein fibres (F-actin and microtubules). They are involved in 
many biological processes essential for living organisms, such as mitosis and meiosis, as well as muscle 
contraction and intracellular transport. 

How individual motors work 

Looking at these systems with the eye of a theoretician, molecular motors reduce to microscopic objects 
(and as such they are subject to Brownian motion) moving uni-directionally along one-dimensional peri- 
odic substrates. The obvious question is: "What are the mechanisms which rectify Brownian motion and 
can explain such a uni-directional motion?" The answer is that there are actually a variety of mechanisms, 
different aspects of which are discussed in a series of recent reviews [237-241] (see also the article by P. 
Hanggi et al. in this volume [242]). 

The essential idea can be summarized with the following rather elementary but illustrative example, 
known as the "flashing" ratchet [243-245]; see Fig. 9. Consider a Brownian particle subject to viscous 
damping which is moving along a track with a periodic, asymmetric potential V(x) (with a sawtooth-like 
shape). Despite the spatial asymmetry of the potential, no preferential direction of motion is possible if in 
an isothermal environment only equilibrium fluctuations act on the particle. This argument goes back to 
Smoluchowski [246] (see also the Feynman lectures [1]) who showed that otherwise one would be able to 
construct a perpetuum mobile and violate the second law of thermodynamics. 
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Fig. 9 Illustration of a "flashing ratchet". Upon switching a sawtooth-like potential on and off, a Brownian particle 
will move uni-directionally along a one-dimensional molecular track. 

Thus, for generating uni-directional motion we obviously need an additional ingredient which drives 
the system out of equilibrium. In a "flashing" ratchet this is achieved by cyclically turning the potential 
on and off. The principle is easily explained. Say that the potential barriers are rather high compared to 
any thermal energy scale. Then, with the potential switched on, the Brownian particle is well localized in 




"Es bleiben uns von den bekannten Naturprozessen noch die organischen Wesen iibrig. ... [Tiere] verbrauchen also eine 
gewisse Quantitat chemischer Spannungskrafte, und erzeugen dafiir Warme und mechanische Krafte." At that time the term 
"Spannungskraft" and "lebendige Kraft" stood for what we nowadays call potential and kinetic energy, respectively. 
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one of the potential minima. Upon turning the potential off, the Brownian particle starts performing a free 
one-dimensional random walk, such that after some time r g the probability distribution P(x) for finding 
it at position x is given by a Gaussian 

P{x) — = exp 

V47rL>T off 

where D is the diffusion coefficient. When the potential is turned on again, the asymmetry of the potential 
acts as a rectifyer on the probability distribution function and causes net transport in one direction; see 
Fig. 9. In this sense Brownian motion acts as a "life force" ("Lebenskraft") for molecular scale engines. 
Note, however, that the energy to drive the transport does not come from the thermal noise but is provided 
when the potential is switched on and off. 

Recently there have been several experimental ramifications of the "flashing ratchet" using colloidal 
particles [247,248]. With the advances in nanosciences and microfluidics there are by now many cleverly 
designed systems, classical as well as quantum, which utilize Brownian motion for generating transport, 
sorting particles and performing other tasks [249]. 

The "flashing ratchet" is, of course, not a realistic model for any kind of actual biological motor. It 
oversimplifies a variety of features of real biological engines. First, molecular motors are not structureless 
particles but proteins with a quite complex architecture and hence a large number of internal degrees of 
freedom. For example, kinesin is a dimer consisting of two globular head domains joined together by alpha- 
helical coiled-coil. Second, the interactions of the motor protein with its fuel ATP and the molecular track 
are complex molecular processes, a full analysis of which would require an atomistic description including 
all the different types of forces between the chemical agents in a watery environment. To go beyond 
generic models making the proof of principle one has to design more elaborate models which account for 
the most important structural and chemical elements. What "most important" means varies from system 
to system. It depends on specific features of the particular engine, e. g. the type of coupling between the 
biochemical cycle and the conformational states of the protein 23 . Depending on the level of coarse-graining 
we may describe the forces acting on the degrees of freedom as systematic or stochastic. Systematic 
forces include viscous drag and "mechanical" forces acting on structural elements of the protein which are 
sometimes visualized as parts of macroscopic engines (springs, levers, joints, etc.). Stochastic forces may 
be classified as Brownian conformational fluctuations and stochastic chemical reactions. Both, of course, 
are just variations of the same theme. Their common cause is the chaotic dance at a molecular scale. The 
"flashing ratchet" and the "Possion stepper" emphasize these two sides of the medal. In the flashing ratchet 
all stochasticity results from Brownian conformational fluctuations in a potential landscape. In contrast, 
for the "Poisson stepper" (a textbook example for a discrete stochastic process), which advances step by 
step at constant rate, stochasticity results from the distribution of waiting times between the steps [250]. 

Traffic jams in the cell 

There are several biological processes where a concerted action of molecular motors is of importance. 
A particulary prominent example is protein sythesis by mRNA translation in eucaryotes, which involves 
unidirectional motion of ribosome complexes along mRNA strands. Theoretical investigations go back to 
the pioneering work by MacDonald et al. [251], who designed a driven lattice gas model which by now 
is known as the totally asymmetric exclusion process (TASEP). In this model a single species of particles 
is hopping unidirectionally and with a uniform rate along a one-dimensional lattice. The only interaction 
between the particles is hard-core repulsion, which prevents more than one particle from occupying the 
same site on the lattice; see Fig. 1 1 . Originally intended as a model for an important biological process, this 
model has by now become one of the paradigms of non-equilibrium physics (for a review see Refs. [252- 
255]). 



x 2 



(26) 



It is impossible to list the multitude of theoretical and experimental work in this area. The reader may consult some of the 
recent reviews [238-241] and books [236] for a more elaborate discussion and additional references. 



24 



E. Frey and K. Kroy: Brownian motion in soft matter and biological physics 



ribosomes' 





Fig. 10 Illustration of protein synthesis: ribosomes 
move from codon to codon along a mRNA template, 
reading off genetic information and thereby generat- 
ing proteins step by step [256] 



Fig. 11 Illustration of the Totally Asymmetric Sim- 
ple Exclusion Process with open boundaries. The en- 
trance and exit rates at the left and right end of the 
one-dimensional lattice are given by a and /3, respec- 
tively. 



Much of what we know about the dynamics and the steady state of the TASEP and driven lattice 
gases [60] in general is the result of theoretical work using a broad variety of tools and exploiting the 
intimate connections between classical and quantum physics, which we eluded to in section 4. It has been 
found that the nature of the non-equilibrium steady state of the TASEP depends sensitively on the boundary 
conditions. For periodic boundary conditions the system reaches a steady state of constant density. Inter- 
estingly, density fluctuations are found to spread faster than diffusively [257]. This can be understood by 
an exact mapping [258] to a growing interface model, whose dynamics in the continuum limit is described 
in terms of the Kardar-Parisi-Zhang equation [259] and its cousin the noisy Burgers equation [260]. In 
contrast to such ring systems, open systems with particle reservoirs at the ends exhibit phase transitions 
upon varying the boundary conditions [26 1 ] . This is genuinely different from thermal equilibrium systems 
where boundary effects usually do not affect the bulk behavior and become negligible if the system is large 
enough. In addition, general theorems do not even allow equilibrium phase transitions in one-dimensional 
systems at finite temperatures (if the interactions are not too long-range) [262]. Yet another difference 
between equilibrium and non-equilibrium processes can be clearly seen on the level of its dynamics. If 
transition rates between microscopic configurations are obeying detailed balance the system is guaranteed 
to evolve into thermal equilibrium. Systems lacking detailed balance may still reach a steady state, but 
at present there are no universal concepts like the Boltzmann-Gibbs ensemble theory for characterizing 
such non-equilibrium steady states. In most instances one has to resort on solving nothing less than its full 
dynamics. It is only recently, that exact (non-local) free energy functionals for driven diffusive systems 
have been derived [263,264]. 

Recently, driven lattice gas models have found their way back into biological physics. MacDonald's 
original ideas have been generalized taking into account the finite extension of ribosomes [265, 266]. Also 
some aspects of intracellular transport, where the molecular motor kinesin moves along microtubules, 
show close resemblance with the TASEP. As an additional feature one has to account for the fact that 
microtubules are embedded in a cellular environment with a reservoir of motors in solution. This allows for 
the motors to attach from the solution to the molecular track and to detach from it to join the reservoir [267, 
268] . Typically kinetic rates are such that these motors walk a significant fraction along the molecular track 
before detaching. It came as quite a surprise that the competition between the non-conservative on/off- 
dynamics and the driven transport along the molecular track leads to a separation between high and low 
density phases [267]. 



adp header will be provided by the publisher 



25 



9 Conclusions and Outlook 

After 100 years many of the central questions surrounding the ubiquitious phenomenon of Brownian mo- 
tion are still unresolved so that it seems still too early for closing the books and drawing final conclusions. 
As soon as we leave the caved paths of equilibrium statistical mechanics traced out a century ago, intrigu- 
ing open problems and promising perspectives abound. Still, we cannot derive the basic assumptions of 
the founding fathers from the underlying microscopic dynamics, and still we have not managed to general- 
ize reliably the powerful theoretical framework we inherited from them to processes far from equilibrium. 
We are only starting to appreciate to what extent the processes of life are shaken by fluctuations, just like 
the pollen grains studied by Brown almost 200 years ago. Life undeniably shows a strong preference for 
the borderline between perfect order and the erratic dance called Brownian motion, which will, for the 
foreseeable future, remain a paradigm of the natural sciences in a noisy world. 
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